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A SIMULATION OF HIGH LATITUDE F-LAYER INSTABILITIES 
IN THE PRESENCE OF MAGNETOSPHERE-IONOSPHERE COUPLING 


I. INTRODUCTION 

Recent experimental ground-based, rocket, and satellite (HILAT, 
DYNAMICS EXPLORER, AURE0L~3) observations in the auroral zone and polar cap 
ionosphere have indicated the existence of both large [Weber et al., 1984; 
Basu et al., 1984; Bythrow et al., 1984; Cerisier et al., 1984; Vickrey et 
al., I960] and small [Hanuise et al . , 1981; 5aker et al., 1983] scale 

density structures and irregularities. Different mechanisms, have been 
proposed to account for high latitude ionospheric irregularities , e.g., 

particle precipitation, plasma instabilities, and neutral fluid turbulence 
[Keskinen and Ossakow, 1983]. Considerable quantitative progress has been 
made in explaining ionospheric structure using plasma interchange 

instabilities, e.g., the Rayleigh-Taylor instability [Balsley et al., 1972; 
Ossakow, 1981] in equatorial spread F, and the E * 3 and current-convective 
instabilities in the high latitude ionosphere [Ossakow and Chaturvedi, 
1979; Keskinen and Ossakow, 1983]. Recently, Weber et al. (1984) and 
Cerisier et al . (1984) have invoked the E x B instability to explain large 
scale density fluctuations in the high latitude ionosphere. 

A shortcoming of most past research on the nonlinear theory of 
the E x 3 instability as it applies to ionospheric structure is that it has 
been restricted to the collisiona-L (or non-inertial ) regime, i.e., vc >> w 
where is the ion-neutral collision frequency and cu is the wave 
frequency. One exception is the recent work of Huba et al . (1985) who 

studied the nonlinear evolution of interchange instabilities in both the 
inertial and non-inertial regimes. However, their work was limited to 

short wavelength turbulence, i.e., kL >> 1 where k is the wavenumber and L 
is the density gradient scale length, which is not applicable to large 
scale ionospheric structures. 

In this letter, we present the first simulation of inertial high 
latitude ionospheric interchange instabilities (e.g., the E x B 
instability) with inclusion of magnetospheric coupling effects. The basic 
conclusions of this study are (1) magnetospheric coupling effects reduce 
the growth rate of the E x B instability, (2) striations produced by the 
inertial E x B instability develop in a different manner than in the non- 
Manuscript approved April 17, 1985. 
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inertial (v., >> to) regime, (3) in configuration space, the striations in 
the inertial regime are more isotropic and spread out resulting in 
irregularities oriented perpendicular to those produced in the non-inertiai 
case . 


II. MODEL 

The physical configuration and assumptions of our model are described 
as follows. We only consider structure in the plane transverse to the 
ambient magnetic field, i.e., the xy plane. The F-layer is initially 
characterized by a 6 to ' density enhancement with a Gaussian profile of 
scale size 12 km in the x-direction and uniform in the y-direction, a 
uniform magnetic field in the z-direction (B_, = 0.5 G), and a background 
electric field in the y-direction (E = .025 V/m). The entire enhance- 
ment E x Bdrifts in the x-direction at a velocity v x = 0.5 km/sec. A 
uniform horizontal magnetosphere is assumed above the F-layer linked by the 
vertical magnetic field lines. The back edge of the F-layer enhancement, 
relative to the drift, is unstable to the E x 3 instability, which drives 
ion Pedersen currents in the F-layer and ion polarization drift currents in 
the magnetosphere due to the perpendicular electric field mapping along the 
geomagnetic field. These currents close along the magnetic field via 
parallel electron currents. Assuming that all ion drifts associated with 
the perpendicular currents are npgligible compared to the E x B drift and 
that E_!_ = -V x 4>, the equations describing this system are 


3n 

— ♦ Vi • nv =0 

(1 ) 

3t -ij_ 

cV 1 4> 

v = v = - — x Z 

(2) 

~ s i. -i-L B 

cnev. 

■If, ’ ' ( bo. ) ‘ °p v x* 

Imj. ■ - - * V 

(3) 

(4) 

0 = / F dz (V -j ) + J M dz (V .J ) 

(5) 



where n is the ion (and electron) density, x e (i)± is the perpendicular 
electron (ion) velocity, and are the F-layer Pedersen ;and the 
magnetospheric polarization drift current densities, respectively, Op is 
the F-layer Pedersen conductivity, V A = B/ ( ) 1 ^ is the Alfven velocity, 
and V| and ^ are the ion-neutral collision frequency and ion gyrofre- 
quency, respectively. Using (2)-(4) in (5) yields the potential equation, 


o - Vj. • (j ♦ C [1- - x z) • JjJIVj.* 

p M 3 1 B 

where the field-line integrated F-layei Pedersen conductivity is 

e d 3 L dz a o’ 


and the field-line integrated magnetospheric inertial capacitance is 

C M*i/ H dZ (’ 2 /V a)' (8 > 

Since the magnetospheric layer is uniform and the ion flow given by (2) is 
incompressible, the magnetosphere remains uniform and the continuity 
equation for this layer may be neglected. The F-layer continuity equation 
for n may be written as a continuity equation for Z , 


3t L p + Vl ’ ^ Vii ] = °‘ 


Therefore, the system is completely described by (2), (6) and (9) in the 
variables l and <j>,- with C H a constant capacitance describing the magnitude 
of magnetosphere-ionosphere coupling. 

The numerical methods used to simulate the model equations are 
described in Zalesak et al. [1982]. The continuity equation (9) is solved 
numerically using the multi-dimensional f lux-corrected techniques of 
Zalesak [1979], while the potential equation (6) is solved with the 
incomplete Cholesky conjugate gradient algorithm c' Hain [1980]. The 
simulations are performed on an 100 x 80 cell grid (x,y) with a cell size 
of 1.0 km x .25 km which is drifting with the enhancement at the E y x B 
velocity. Periodic boundary conditions are assumed in the y-directicn, and 
the grid is initialized with a random density fluctuation. 



If Ep has scale length L~ = (l /E p }3Ep/9x, the non-inertial (C M = 0) 
growth rate for a wave with wave vector k parallel to y is = cEy/BL in 
the regime kL >> 1 . For the parameters of our simulation, the maximum (non- 
inertial) growth rate is Yq =* 0.05 sec 1 . The e-folding distance 
of Ej_ parallel to B for this wave is about Xj - (o | / oj_ ) 1 7 ^kj_ 1 [Volk and 
Haerendel, 1971], implying that a wave with k_L _1 =» 1 km has an e-folding 
distance of 10-^ km parallel to the field and that Ej_ maps well up into the 
magnetosphere where the e-folding distance is much greater due to the rapid 
decrease in a . The Alfven speed is V fl - 10^ km/sec so that electric 

P H 

fields of a wave with a growth time of 20 sec map upwards on the order of 

4 1 7 

10 km into the magnetosphere, and C M has a value of roughly 10 J cm 

(or - 10 farad). 

The effect of ion inertia on the E x 3 instability has previously been 
studied by Ossakow et al . [1978] in the linear regime. Although their 
results were for Pedersen and oolarization drifts at the same altitude, the 
linear growth rates are effectively the same for our model equations. A 

r.on-zero inertia implies the existence of an inertial relaxation rate, 

Y i = Ep/C M . The growth rate in the presence of inertia for the short 
wavelength approximation has two regions of interest, the non-inertial 
regime 

Y * Yq for 4 Yq << Y. (10) 

and the inertial regime 

Y * (YqY.) 172 for 4 y q >> Y.. (11) 

Assuming a E for the F-layer of 3 x 10 10 cm/sec (- 0.1 mho), a typical 
^ _ 1 

value of Y i for the simulation is .008 sec , which is in the inertial 

regime and which reduces the growth rate of the instability from .05 sec 1 

to .02 sec . For our simulation, two cases were run: a non-inertial case 

= 1.00 sec’ 1 , and an inertial case for which Yj = 0.01 sec . 


for which Y.- 


III. RESULTS 


The results of the two simulations are shown in Figures 1 and 2. The 
four panels in each simulation show approximately equal times relative to 
the linear growth time of the instability. Figure 1 shows the results for 
the non-inertia] E x g instability, while Figure 2 shows the results for 
the inertial case. It is clear that the results of the two cases are very 
different at later times. 

The behavior of the plasma in Figure 1 is typical of that observed in 
previous E x B instability simulations [Keskinen and Ossakow, 1982]. In 
panel 2, at 104 seconds, we see that early in the nonlinear stage a set of 
"fingers" has clearly formed. The high density fingers grow outward into 
the low density plasma while the low density fingers penetrate into the 
high density cloud. Subsequent nonlinear evolution involves the continued 
elongation of these fingers with very little apparent change in the size of 
the structures perpendicular to their long dimension as seen in panels 3 
and 4. The original density enhancement has effectively been sliced into a 
group of sheets parallel to the initial density gradient. 

The behavior of plasma and structure formation in Figure 2, the 
inertial case, is very different. During the time period between panel 1 
and panel 2, before 300 seconds, the growth of the instability closely 
resembles that in Figure 1 . There is nonlinear development of long narrow 
high and low density fingers which move in opposite directions. However, 
after 300 seconds, the behavior changes radically. In panels 3 and 4, we 
see that the fingers form mushroom-like heads and tend to thicken. No 
longer are there long thin interpenetrating fingers, rather they are fat 
interpenetrating blocs. Any narrow fingers which begin to form quickly go 
to a mushroom shape and then spread cut. In a number of simulations we 
have noted a tendency for the structure in the > direction to undergo an 
inverse cascade to the longest mode which will fit in the system. This 
feature can clearly be seen in panel 4, where the structured state 
throughout most of the simulation region shows two blobs; one of high 
density, the cthe' 1 of low density. 

The non-inertial simulation, shown in Figure 1 , is characterized by 
striations with a large value of ky, which grow in the x-direction without 
hindrance, completely destroying the initial orientation of the density 
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enhancement in the y-direction. The inertial simulation, shown in Figure 
2, is character ized initially by striations with a large value of ky, which 
are hindered ir: their growth in the x-direction. The plasma at the leading 
edge of the striation finger is swept to either side and around the 
striation forming the characteristic mushroom shape. Further evolution 

involves an inverse cascade in ky to the minimum value allowed on the 
numerical grid while simultaneously, the typical value of k v increases. 

A 

Consequently, individual enhancements tend to an orientation similar to the 
original x-directed orientation. Furthermore, the individual enhancements 
appear to be stabilized to further S x g structuring owing to velcrity 
shear along their apparently unstable faces (Perkins and Doles, 1975; Huba 
et al., 1983). 

Interestingly, Rino et al. (1979), Livingston et al. (1982), and Rit 
and Vickrey (1982) have reported sheet-like structures in the night side 
auroral region which are aligned perpendicular to the large scale F region 
ionospheric density gradient. This alignment perpendicular to the density 
gradient is not easily explained by the traditional instability theory. 

However, the alignment is similar to the final alignment observed in the 

inertial simulations. 

IV. CONCLUDING REMARKS 

Several conclusions can be drawn from our simulations. First, 
instability-generated electric fields in the high-latitude F-layer may map 
well up into the magnetosphere, resulting in a much reduced linear growth 

rate due to the effectively increased ion inertia. Second, the nonlinear 

development of the instability is fundamentally different in the presence 
of this coupling; the striations produced are spread and retarded by ion- 
inertial effects resulting ir. more isotropic irregular ities than in the 
non-inertial case. The inertial effects may even tend to stabilize the 
final nonlinear state by producing a velocity shear across normally 
unstable gradients. The simulation results demonstrate some interesting 
features which may have been observed in the high latitude nighttime 
ionosphere. 
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A plot of density for four times for the ease of 1.00 

sec -1 (the non-inertial case). The shaded region represents those 
areas whose density is greater than 2.5 times the background 
density. The simulation grid is periodic in the direction of E 
and is moving with the E x B velocity. 
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A plot of density for four times for the ease of Y^= 0.01 see 1 
(the inertial case). The shaded region and the simulation grid 
are the same as in rig. 1 . The times shown in the two figures are 
roughly equal multiples of the linear growth time for tne 
instaoility in the two oases. 
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1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
0 1 C Y ATTN W. HALL 
01 CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP 
12340 SANTA MONICA BLVD. 

LOS ANGELES, CA S0025 

01 CY ATTN E.C. FIELD, JR. 


RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA 90406 
01 CY ATTN CULLEN CRAIN 
0 1 C Y ATTN ED BEDROZIAN 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY , MA 01 776. 

01CY ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
330 WEST 42nd STREET 
NEW YORK , NY 1 0036 

01CY ATTN VINCE TRAPANI 
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SCIENCE APPLICATIONS, INC. 

1 1 50 PROSPECT PLAZA 
LA JOLLA, CA 92037 

01CY ATTN LEWIS M. LINSON 
01 Cl ATTN DANIEL A. HAMLIN 
01CY ATTN E. FRIEMAN 
01 CY ATTN E.A. STRAKER 
01CY ATTN CURTIS A. SMITH 

SCIENCE APPLICATIONS, INC 
1710 COODRIDCE DR. 

MCLEAN, VA 22102 
0 1 C Y J. COCKAYNE 
01CY E. HYMAN 

SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 9^025 


01 CY 

ATTN 

J CASPER 

01CY 

A~In 

DONALD NEILSON 

01 CY 

at: 

aN BURNS 

01CY 

A 1 ' 1 

SMITH 

01CY 

AI 

~. TSUNODA 

0 1 C Y 

A t T ! 

AVID A. JOHNSON 

01CY 

/ i'T-1 

ALTER G. CHESNUT 

0 ICY 

ATTN 

CHARLES L. RINO 

01 CY 

ATTN 

WALTER JAYE 

01CY 

ATTN 

J. VICKREY 

01 CY 

ATTN 

RAY L. LEADABRAND 

01CY 

ATTN 

G. CARPENTER 

01 CY 

ATTN 

G. PRICE 

01CY 

ATTN 

R. LIVINGSTON 

01 CY 

ATTN 

V. GONZALES 

0 ICY 

ATTN 

D. MCDANIEL 


VISIDYNE 

SOUTH BEDFORD STREET 
BURLINGTON , MASS 01803 
01 CY ATTN W. REIDY 
01CY ATTN J. CARPENTER 
01 CY ATTN C. HUMPHREY 

UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA 15213 

01CY ATTN: N. ZABUSKY 

DIRECTOR OF RESEA RCH 
U.S. NAVAL ACADEMY 
ANNAPOLIS, MD 21402 
02CY 


TECHNOLOGY INTERNATIONAL CORP 
75 WIGGINS AVENUE 
BEDFORD, MA 01730 

01CY ATTN W P. BOQUIST 

TOYON RESEARCH CO. 

P.O. Box 6890 
SANTA BARBARA, CA 93111 
0 1 C Y ATTN JOHN ISE, JR. 

01 CY ATTN JOEL GARBARINO 

TRW DEFENSE A SPACE SYS GROUP 

ONE SPACE PARK 

REDONDO BEACH, CA 90278 


01 CY 

ATTN 

R. 

K. °LEBUCH 

01CY 

ATTN 

S. 

ALTSCHULER 

0 1 C Y 

ATTN 

D . 

DEE 

0 1 C Y 

ATTN 

D/ 

STOCKWELL 
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NAVA.. RESEARCH LABORATORY 
WASHINGTON, D.C. 20375 
Code 4100 

Dr. H. GURSKY - CODE 4100 
Dr. P. GOODMAN - CODE 4180 

A.F. GEOPHYSICS LABORATORY 
L.G. HANSCOH FIELD 
BEDFORD, MA 01731 
DR. T. ELKINS 
DR. W. SWIDER 
MRS. R. SAGALYN 
DR. J.K. FORBES 
OR. T.J . KENESHEA 
OR. W. BURKE 
DR. H. CARLSON 
DR. J. JASPERS 
Dr. F.J. RICH 
DR. N. MAYNARD 

BOSTON UNIVERSITY 
DEPARTMENT OF ASTRONOMY 
BOSTON, MA 02215 
DR. J. AARONS 

CORNELL UNIVERSITY 
ITHACA, NY 14850 
DR. W.E. SWARTZ 
DR. D. FARLEY 
DR. M. KELLEY 

HARVARD UNIVERSITY 
HARVARD SQUARE 
CAMBRIDGE, MA 02138 

dr. m.b. Mcelroy 

DR. P. LINDZEN 

INST1TUT £ FOR DEFENSE ANALYSIS 
400 ARM' /NAVY DRIVE 
ARLINGTON, VA 22202 
DR. E. BAUER 

MASSACHUSETTS INSTITUTE OF 
TECHNOLOGY 
PLASMA FUSION CENTER 
LIBRARY, NU16-262 
CAMBRIDGE, MA 02139 


NASA 

GODDARD SPACE FLIGHT CENTER 
GREEH8ELT, MD 20771 
DR. K. MAEDA 
DR. S. CURTIS 
DR. N. DUBIN 

COMMANDER 

NAVAL AIR SYSTEMS COMMAND 
DEPARTMENT OF THE NAVY 
WASHINGTON, D.C. 20360 
DR. T. CZUBA 

COMMANDER 

NAVAL OCEAN SYSTEMS CENTER 
SAN DIEGO, CA 92152 

MR. R. ROSE - CODE 5321 

N0AA 

DIRECTOR OF SPACE AND 
ENVIRONMENTAL LABORATORY 
BOULDER, CO 80302 
DR. A. GLENN JEAN 
DR. G.W. ADAMS 
DR. D.N. ANDERSON 
OR. K. DAVIES 
DR. R.F. DONNELLY 

OFFICE OF NAVAL RESEARCH 
800 NORTH QUINCY STREET 
^kuNGTON, VA 22217 
OR. G. JOINER 

PENNSYLVANIA STATE UNIVERSITY 
UNIVERSITY PARK, PA 16802 
DR. J.S. NIS8ET 
OR. P.R. R0HR3AUGH 
DR. L.A. CARPENTER 
DR. M. LEE 
DR. R. DIVANY 
DR. P. BENNETT 
DR. F. KLEVANS 

SCIENCE APPLICATIONS, INC. 
1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 
DR. O.A. HAMLIN 
DR. E. FRIEMAN 


STANFORD UNIVERSITY 
STANFORD, CA 94305 
DR. P.M. BANKS 

U.S. ARMY ABERDEEN RESEARCH 
AND DEVELOPMENT CENTER 
BALLISTIC RESEARCH LABORATORY 
ABERDEEN, ND 

DR. J. HEIMERL 

GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIRBANKS, AK 99701 
DR. L.E. LEE 

UNIVERSITY OF CALIFORNIA, 
BERKELEY 

BERKELEY, CA 94720 
DR. M. HUDSON 

UNIVERSITY OF CALIFORNIA 
LOS ALAMOS SCIENTIFIC LABORATORY 
J-10, MS-664 
LOS ALAMOS, NM 37545 
DR. M. PONGRATZ 
DR. 0. SIMONS 
DR. G. 8ARASCH 
DR. L. DUNCAN 
DR. P. BERNHARDT 
DR. S.P. GARY 

UNIVERSITY OF MARYLAND 
COLLEGE PARK, MD 20740 
DR. K. PAPADOPOULOS 
DR. E. OTT 

JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 
DR. R. GREENWALC 
DR. C. MENG 

UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA 15213 
DR. N. ZABUSKY 
DR. M. BIONDI 
DR. E. OVERMAN 


UNIVERSITY OF TEXAS 
AT DALLAS 

CENTER FOR RESEARCH SCIENCES 
P.O. BOX 688 
RICHARDSON, TX 75080 
DR. R. HEELIS 
DR. U. HANSON 
DR. J.P. McCLURE 

UTAH STATE UNIVERSITY 
4TH AND 8TH STREETS 
LOGAN, UTAH 84322 
DR. R. HARRIS 
DR. K. BAKER 
DR. R. SCHUNK 
DR. J. ST. -MAURICE 

PHYSICAL RESEARCH LABORATORY 
PLASMA PHYSICS PROGRAMME 
AHMEDABAD 380 009 
INDIA 

P.J. PATHAK, LIBRARIAN 

LABORATORY FOR PLASMA AND 
FUSION ENERGY STUDIES 
UNIVERSITY OF MARYLAND 
COLLEGE PARK, MD 20742 
JHAN VARYAN HELLMAN, 

REF LIBRARIAN 
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to replicate rrqolith samples that contain agglutinates, shock-damaged particles, 
and surface-correlated volatiles. How much expense and effort are justified in 
simulation remains to be determined. 

Although some lunar samples will have to be sacrificed for utilization experiments, 
this sample lass must be minimized. Experience gained from the last sixteen years 
of lunar sample research will contribute considerably in this effort For example, 
consortium studies have been exceptionally productive in bringing diverse disciplines 
together for studying complex lunar samples; in utilization experiments, a benefit 
of consortium studies would be the extended life of the sample collection. Sample 
life may be prolonged by sequential experimentation; for example, experiments on 
delicate features such as surface-correlated volatiles might precede melting of the 
same sample for experiments on metal and oxygen production. Innovative research 
programs must be developed, and will require close liaison between experimenters 
with different goals, improved techniques are necessary to obtain critical information 
from samples much smaller than these normally used in terrestrial studies. 

APPLICATIONS WORKSHOP 

Communication and cooperation between researchers is vital An Applications 
workshop is being planned, probably to be convened in 1986, that will help to 
explore research opportunities, to develop an understanding of what can and should 
be simulated, and to open communications between research groups. Major goals 
of this workshop would be to (1) encourage experimenters to form consortia, (2) 
expose utilization experimenters to the available scientific data base on lunar samples, 
and (3) expose the present lunar sample research community to the data needs 
of utilization experimenters. 

This workshop will probably include a general session on the data available 
for lunar rocks, minerals, soils, and volatile elements and will cover the prospects 
for producing useful simulants. Topical sessions will discuss fuel extraction (oxygen, 
hydrogen aluminum, and calcium), metal extraction methods (chemical electrical, 
thermal and mechanical), and metal fabrication and the formation of structural 
materials (concrete, glasses, ceramics, and sintered products). A session on mining 
and bulk materials processing will encompass regional geology, strip and subsurface 
techniques, ore processing and the estimation of probable ores. A session on lunar 
environmental control and protection will cover the topics of dust waste, and 
atmosphere control. A special session must be reserved for fresh ideas about other 
innovative products and materials for the space environment 
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RELATIONSHIP 
TO MISSIONS 


Lunar sample studies support all NASA planetary missions in a general way 
by elucidating planetary processes such as volcanism, impact cratering, and planetary 
differentiatioa Sample studies relate directly to the Lunar Geoscience Observer, and 
possibly to other planned missions as well. 

LUNAR GEOSCIENCE OBSERVER (LGO) 

The Lunar Geoscience Observer might follow the Mars Observer in the Planetary 
Observer series Lunar research has provided well-defined scientific goals for this 
mission, i.e., to resolve the problems outlined here in the section entitled Outstanding 
Lunar Science Problems. Scientists familiar with lunar samples are actively involved 
in planning the missioa LAPST intends this involvement to continue into flight teams. 
Lunar sample results are and will continue to be the primary source of "ground 
truth" data that will ensure the fullest use of the geochemical and mineralogical 
global maps to be provided by the LGO missioa Because the upper lunar crust 
is extensively mixed by impacts, continued sample study will refine and extend the 
spectrum of lunar primary igneous rock types and is critical to our understanding 
of regional geochemical compositions, which will undoubtedly be interpreted as 
mixtures of pristine rock types. The discovery of pristine rock types could greatly 
alter our perception of the geologic and petrologic significance of geochemical 
variations that will be mapped by LGO. New laboratory work is needed on the 
spectral characteristics of a variety of lunar rock types in order to interpret the 
global maps of lunar spectral reflectance that LGO will obtain. Earth-based telescopic 
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spectra have shown us which lunar rock types are potentially relevant to the 
understanding of lunar spectral data. The larger the spectral data base on lunar 
samples that is available prior to the LGO mission, the better we will be able to 
comprehend the abundant and rapidly-returning spectral data from the spacecraft 
during mission operations The result: a rvrw understanding of and perspective toward 
the Moon's composition and evolution that will most likely feed back to more specific 
sample studies 

SAMPLE RETURN FROM THE MOON AND OTHER BODIES 

Sample return missions to various large bodies in the solar system are highly 
desirable because some fundamental questions can only be answered by making 
measurements on samples These bodies include Mars, the Moon, and Venus However, 
the rationale for sample returns needs to be addressed. Below are some of the 
key points that demonstrate the advantages of sample return over in situ 
measurements 



Oblique view of the Aristarchus Plateau (PLj and environs The fresh. Copemican-age crater Aristarchus (A) is 
about 42 km in diameter, remote-sensing data indicated the presence of dinopynaene- bearing high Th /—IS 
ppm) rocks —possibly KREEP-nch. plutonic evolved rocks such as granite or quartz monzodiorite The Pla- 
teau itself displays varied geological units including the Imbnum Basin Alpes Formation (an ejecta fades), 
light plains (P) for which Th contents suggest KREEP lava flows, and extremely red. KREEP-rich pyroclastic 
dark mantle deposits Schmid's Valley (S) is a complex, nested sinuous rille. with a morphology that implies a 
protracted volcanic eruption history Mare basalts (M) on the horizon span an age from about SO by to possi- 
bly Jess than 10 by. A future sample return from this area of the Moon has the potential to collect a wide var- 
iety of rock types within a relatively limited geographic region 
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Sample Return Rationale 

Laboratory Equipment Will Always Have Better Resolution and Precision Than 
Flight Instruments. Over the past ten years, laboratory instruments have 
undergone major improvements in their ability to analyze increasingly smaller samples 
and to provide significantly greater resolution of physical and chemical differences 
and improved precision for isotopic analysis. Because of weight and power limitations, 
no proposed flight instrument can ever be as good as the best laboratory instrument 
in its ability to analyze small samples or to detect small differences. This ability 
is critical to understanding the samples, so laboratory instruments will always provide 
better data and better understanding of planetary samples. For many types of planetary 
samples this superior resolution and sensitivity is not only desirable but is absolutely 
vital to understanding the material. For instance, only by using instruments of very 
high resolution that can discriminate among individual mineral and glass phases 
is it possible to understand a planetary regolith. Very precise data are needed for 
reliable isotope age determinations and isotopic ratios. It is not likely that some 
kinds of analyses will ever be accomplished to the required precision by remote 
instruments. Such age and isotopic ratio determinations are critical to the 
understanding of the geological evolution of a planet and can be used to calibrate 
other time scales such as those based on crater densities. 

Returned Sar. pies Become Resources That j re Accessible in the Future for Rapidly 
Improving Analytical Technology. Over a period of years, analytical instruments 
undergo considerable improvements and entirely new instruments are developed 
If planetary samples are brought back to terrestrial laboratories, improved or new 
instruments can be applied to them as they become available. The samples can 
then be periodically mined for new data. In contrast, flight instruments 'freeze’’ the 
state of the technology some time before launch and cannot be improved Only 
an entirely new mission can take advantage of improvements in instrument technology 
that may occur from year to year. The technology of mass spectrometry could not 
support Nd-Sm dating at the time of the Apollo missions, but several years later 
this technique provided key data on returned lunar samples. Third-generation ion 
probes and PIXIE analyses now allow trace element determinations on individual 
mineral and glass grains; such techniques were not available during the Apollo 
missions. 

Unforeseen Key Discoveries Can Lead to New Experimental Design. A package 

of flight instruments must be chosen and designed well before the mission begins. 
The choice of the instrument mix and the design of each instrument is entirely 
dependent on the best guess of what the targeted samples may be like. If the samples 
are not entirety as expected, the instruments may not work at all, may give the 
wrong kind of data, or may miss the critical data. An example is an instrument 
package designed to sample fine particulate material that instead encoun^rr.1 coarse- 
grained material. If the chemistry of the samples is not as expecteu, the proper 
analyzing system may not be included in the package. Flight instruments inevitably 
lack flexibility and can only do what they are designed to do based on a model 
for the characteristics of the samples to be analyzed. If the sample characteristics 
were known perfectly before the flight, the analyses would not be necessary. He w.--, 


the supposed characteristics are likely to be wrong in some respects and, therefore, 
the package will not be entirely suitable. Sample return and laboratory instruments 
allow complete flexibility of instrument types and sample preparation, and even 
have the ability to handle totally unexpected situations 

There are No Weight/Power Considerations on Laboratory Instruments Many 

laboratory techniques require massive or power-hungry devices. Examples include 
high-resolution mass spectrometers, ion probes, ferromagnetic and nuclear magnetic 
resonance devices, and synchrotron radiation. A recent example is the Brookhaven 
synchrotron XRF probe for lunar and meteorite samples. None of these techniques 
are easily adaptable to flight instruments because of the power and weight limitations. 

The Variety and Complexify of Laboratory Instrumentation for Sample Studies 
are Unlimited. Laboratories contain a variety of very complex instruments. Literally 

hundreds of different types of instruments and techniques have been applied to 
lunar samples and have provided a complex variety of data types, all of which 
have contributed toward understanding the samples. In contrast, only a few different 
instruments will ever be practical for a flight package. Even if these instruments 
are equal to their laboratory equivalents, the limited number of instruments will 
limit the kinds of data that can be acquired. If we had been limited to only ive 
or ten instruments in our analyses of lunar samples, it is doubtful that we would 
have made much progress in understanding these complex samples, in addition, 
sterilization requirements can limit the quality and complexity of flight instruments. 

OTHER POSSIBILITIES 

A central part of NASA's program is the construction of the Space Statioa 
It is difficult to predict how the Space Station will affect lunar science, but we envision 
several possibilities. One is improved telescopic observations of the Moon from an 
observatory installed on or near the Space Statioa The absence of atmosphere 
would allow observations of the Moon in a wider range of wavelengths than is 
possible from Earth's surface, thus providing more chemical and mineralogical 
information. Such an observatory would supplement the LGO. 

Experiments in microgravity and high vacuum could have significance for 
understanding lunar rock evolution. Such experiments have not yet been designed, 
but we expect them to be proposed when material research laboratories are in 
place on the Space Statioa The Space Station will be the best place to conduct 
many lunar-sample utilization studies, for many lunar materials will ultimately be 
used in orbit for space manufacturing and construe tioa Again, there is interaction 
progression between sample studies and sample utilizatioa 


OUTLOOK 

FOR CURATORIAL OPERATIONS 


It is vital that the curatorial operations continue to support the Planetary Materials 
and Geochemistry Program (PMGP), at least at the present level of activity. In 1981, 
the curatorial effort was significantly reduced and has since remained at a constant 
level of effort, with budget growth sufficient to cover inflationary costs only. The 
present activity is efficiently supplying the needs of the community and preserving 
the integrity of the samples. 

The active approach of exposing more surfaces and samples for study by the 
planetary materials community should be continued Cutting of breccias is stimulating 
significant research activities Several research groups have also taken advantage 
of the chance to pick significant samples from the supplies of 1 -4 mm soil fragments. 
However, detailed cataloging would permit more effective use of these valuable 
samples. Systematic thin-sectioning and descriptions of regolith breccias have focused 
the research of several groups on these samples. LAPST will continue to encourage 
the Curator to actively process the collectioa 

LAPST will evaluate the need to open specific lunar cores for investigatioa 
The dissection of lunar cores was discontinued in 1981 with approximately one- 
third of the cores unexamined Several lunar investigation teams have expressed 
interest in studying new cores. LAPST will evaluate both the significance of the 
potential science and the impact on curatorial operations. 

LAPST will continue to review the status of the cataloging and documentation 
of the lunar sample collectioa The present set of documentation is uneven in 
thoroughness and quality Excellent, thorough reviews in the form of catalogs are 
available for the Apollo 15 and 16 collections LAPST will systematically review 


the other sample documentation (guidebooks, catalogs, core descriptions, etc.) and 
identify significant shortcomings. A prioritized plan will be proposed to direct such 
future curatorial efforts. 

LAPST and the Curator will remain ready to take advantage of opportunities 
to develop new sources of extraterrestrial materials. For example, recently returned 
Solar Max spacecraft parts have provided samples of captured cosmic dust The 
Long Duration Exposure Facility (LDEF) will soon be available as another source. 
In the new era of replaceable, repairable spacecraft, LAPST expects new opportunities 
and will encourage use of the curatorial faclities to recover cosmic dust samples. 

A significant challenge to the Curator and LAPST is the increasing interest in 
studies of potential utilization of lunar samples and other extraterrestrial materials. 
The Planetary Materials Branch is the sole source of lunar materials, and the curatorial 
staff will identify samples that are suitable for use in such engineering and applied 
studies. The curatorial staff will be a source of information on suitable simulants 
for engineering and applied studies. Use of simulants can reduce the demand for 
lunar material and should be used to pre-test experiments in order to optimize 
return from studies of actual lunar materials. LAPST recognizes the significant costs 
of producing simulants but also recognizes the need to have studies done on realistic 
lunar simulants. Whether major stocks of simulated lunar material should be produced 
by the Curator is a question that cannot be answered until more is known about 
the nature, volume, and cost of responding to requests for such materials. 
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